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ABSTRACT 

We present high-resolution spectroscopy of the neutron star/low-mass X-ray binaries (LMXBs) 
4U 1850—087 and 4U 0513—40 as part of our continuing study of known and candidate ultracompact 
binaries. The LMXB 4U 1850—087 is one of four systems in which we had previously inferred an 
unusual Ne/O ratio in the absorption along the line of sight, most likely from material local to the 
binaries. However, our recent Chandra X-ray Observatory LETGS spectrum of 4U 1850—087 finds 
a Ne/O ratio by number of 0.22±0.05, smaller than previously measured and consistent with the 
expected interstellar value. We propose that variations in the Ne/O ratio due to source variability, as 
previously observed in these sources, can explain the difference between the low- and high-resolution 
spectral results for 4U 1850-087. Our XMM-Newton RGS observation of 4U 0513-40 also shows no 
unusual abundance ratios in the absorption along the line of sight. We also present spectral results 
from a third candidate ultracompact binary, 4U 1822—000, whose spectrum is well fit by an absorbed 
power-law + blackbody model with absorption consistent with the expected interstellar value. Finally, 
we present the non-detection of a fourth candidate ultracompact binary, 4U 1905+000, with an upper 
limit on the source luminosity of < 1 x 10 32 erg s _1 . Using archival data, we show that the source 
has entered an extended quiescent state. 

Subject headings: binaries: close — stars: individual (4U 0513-40, 4U 1822-000, 4U 1850-087, 
4U 1905+000) — X-rays: binaries 



1. INTRODUCTION 

Low-mass X-ray binaries (LMXBs) consist of a neu- 
tron star (NS) or black hole (BH) accreting material 
from a <1 M Q donor star. One particularly interest- 
ing subset of the LMXB population has orbital peri- 
ods <80 min, the predicted minimum period for bina- 
ries with hydrogen-rich main sequence donors. Such sys- 
tems, commonly called ultracompact binari es, require 
hydrogen-deficient or degenerate companions lloss et al.l 
119781 iRappaport fc Josslll984l iNelson et al.lll986[) . Out 
of the 53 LMXBs with orbital period measurements, 
eight systems have ultracompact periods. The growing 
number of known and candidate ultracompact binaries 
suggests that these systems make up a significant frac- 
tion of the LMXB population. We have undertaken X- 
ray spectral observations of one known and two candi- 
date ultracompact binaries to better understand these 
unusual systems. 

The LMXB 4U 1850-087 is an X-ray burster 
located in the globular cluster NGC 6712. It 
detected by all of the major X-ray satellites 
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19801: IWarwick et alJll981i iPriedhorskv fc TerreL.. 
Wood et alJ 119841: IWarwick et alJ 11988b iKitamoto et all 
1992HChristian fc SwanldM 997H.Iuett et a.Ul2001^ . Dur- 
ing the 1995 ASCA observation, the 0.5-10 keV X-ray 
flux was 3.4 x 10~ 10 erg cm~ 2 s _1 l| Juett et al J 12001^ . 
which is representative of the average source flux. A low- 
amplitude 20.6-min periodicity, attributed to the orbital 
period, was reported in the HST data of the ultraviolet 
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was 

since Uhuru and shows a roughly factor of 
variation in flux llForman et alJ IT978; 



counterpart of 4U 1850-087 l|Homer et al.lll996|) . While 
this periodicity has yet to be confirmed, the X-ray to 
optical luminosity ratio (Lx/L opt ) is consistent with an 
ultracompact binary based on the relationship between 
optical magni tude, X-ray luminosity, and orbita l period 
developed bv Ivan Paradiis fc McClintockl lfl99l . Addi- 
tionally, the X-ray luminosity of 4U 1850—087 is con- 
sistent with the expected mass-transfer rate for a sys- 
tem with an orbital period of 21 min assuming a WD 
donor and gravitat ional radiation driven mass transfer 
ijHomer et al.lll996|) . 

Spectral evidence als o supports an ultr acompact na- 
ture for 4U 1850-087. iSidoli et all 1)200 lj) analyzed the 
BeppoSAX spectra of the majority of the bright LMXBs 
located in globular clusters. They found that the best 
fit spectral parameters seemed to fall into two groups, 
with one group consisting of 4U 1820-30, 4U 1850-087, 
and 4U 0513—40, i.e., the known or candidate ultracom- 
pact binaries. Additionally, we found that the ASCA 
spectrum of 4U 1850—087 revealed an unusual feature at 
0.7 keV also seen in three other candidate ultracompact 
binaries l) Juett et al.1 12001)) . We attributed this feature 
to an unusual Ne/O ratio in the interstellar absorption 
edges, presumably influenced by material local to the 
binary. However, a recent observation of 4U 1850—087 
using XMM -Newton found no evidence for an unusual 
Ne/O ratio (jSidoli et al.ll2004)L 

The LMXB 4U 0513-40 is an X-ray burster in the 
globular cluster NGC 1851. The binary was first dis- 
covered by the MIT OSO-7 satellite and has been stud- 
ied by all of the major X-ray instrument s that followed 
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ll999t ISidoli et al.l|20T)l . The flux from 4U 0513-40 
shows a factor of ten variability from observation to 
observation. During the 2000 BeppoSAX observation, 
4U 0513-40 ha d a 0.1-100 keV X-ray luminosity of 
5 x 10 36 er e; s" 1 (j Sicioli et aIll2Q01l) . an average flux for 
the source. iDeutsch et all l|200Cft identified 4U 0513-40 
as a candidate ultracompact binary based on its Lx/L opt 
ratio. In addition, its broadband optical, ultraviolet, 
and X-ray sp ectrum resembles that of other ultrac om- 
pact LMXBs ijDeutsch et al.ll2000t ISidoli et al.ll2001(l . 

We identified the Galactic field LMXB 4U 1822-000 
as an ultracompact binary candidate based on its 
Lx/L pt ratio. The source was discovered by Uhuru 
and has be en observed by all of the major X- 
ray missions llGiacconi et all [1973 iMarkert et al J 119711 
I Warwick et al.ll98lUWood et al.ll984HSmale fe Lochnerl 
119921 IGottwald et alJ 119951 IChristian fc Swankl 11997^ 
The source 4U 1822-000 has been roughly constant in 
flux over the course of X-ray astronomy, with a 0.5- 
20 keV unabsorbed flux of 1 x 10~ 9 erg cm~ 2 s _1 
ijChristian fc Swanklll997|) . It has a faint optical coun- 
terpart ( V=22) similar to other cand idate ultracompact 
binaries l)Chevalier fc Ilovaiskvlll985j) . 

In the published literature, 4U 1905+000 was found 
to be persistently bright, with a 2-10 keV flux 
of 3 x 10 -10 erg cm" 2 s" 1 , and showed Type-I 
X-ray bursts (ILewin et ail 1197a ISeward et all 119761 
Forman et all Il978t iReid et alJ Il980t IChevalier et all 
19851 IChevalier fc Tlovalskvlll 9901: IGottwald et al.Nl 9951 
Christian & Swank 119971) . These observations were 
performed in the 1970s and 1980s. The source was 
observed by more recent X-ray satellites, but the 
results of these observations have never been pub- 
lished. From the ana l ysis o f a radius expansion burst, 
IChevalier fc Ilovaiskvl l)1990() estimated the distance to 
4 U 1905+000 at 8±1 kpc . The optical counterpart found 
bv lChevalier fc Ilovaiskvl <| 19851) is a V=20.5, ultraviolet- 
excess star. 

In this paper, we present the results from X-ray spec- 
tral observations of 4U 1850-087, 4U 0513-40, and 
4U 1822-000. The non-detection of 4U 1905+000 is also 
presented. The details of the observations and data re- 
duction techniques are given in §21 an d the analysis is 
presented in 331 In we discuss the implications of 
these results and summarize the case for neon-rich donors 
in ultracompact binaries. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Chandra Observation of 4U 1850-087 

We observed 4U 1850-087 with the Chandra X- 
ray Observatory on 2002 August 16 for 50 ks using 
the Low Energy Transmission Grating Spectrometer 
(LETGS) and the Advanced CCD Imaging Spectrom- 
eter ("ACTS: fWeisskonf et, abll200l . The LETGS spec- 
tra are imaged by ACIS, an array of six CCD detectors. 
The LETGS/ ACIS combination provides both an undis- 
persed (zeroth order) image and dispersed spectra from 
the grating with a first order wavelength range of 1.4- 
63 A (0.2-8.9 keV) and a spectral resolution of AA = 
0.05 A (i?=A/AA=28-1260). The various orders overlap 
and are sorted using the intrinsic energy resolution of 
the ACIS CCDs. The observation of 4U 1850-087 used 
a Y-offset of +1.'5 in order to place the O-K absorption 



edge on the back-side illuminated S3 CCD, which has 
suffered less degradation than the front-side illuminated 
CCDs. The detected zeroth and first order count rates 
for 4U 1850—087 were 1.2 and 3.0 cts s _1 , respectively. 

The "level 1" event file was processed using the CIAO 
v3.0 data analysis package 2 . The standard CIAO spec- 
tral reduction procedure was performed. We filtered 
the event file retaining those events tagged as afterglow 
events by the acis_detect_af terglow tool. Since order- 
sorting of grating spectra provides efficient rejection of 
background events, the afterglow detection tool is not 
necessary to detect cosmic ray afterglow events. No fea- 
tures were found that might be attributable to afterglow 
events. For bright sources, pile up ca n be a problem for 
CCD detectors (see e.g.. iDavisl 12003). The zeroth order 
Chandra spectrum of 4U 1850—087 was heavily affected 
by pileup and was not used in this analysis. The first or- 
der spectrum suffered from minimal pileup in the 6-8 A 
range (1.5-2.1 keV) which has the effect of making the 
instrumental iridium edge feature more pronounced. De- 
tector response files (ARFs and RMFs) for the plus and 
minus first order spectra were created using the standard 
CIAO tools. 

Calibration studies have determined that a layer of 
contaminant is building on the optical blocking filter 
of ACIS, absorbing incoming X-rays. The contami- 
nant gradually reduces the detector quantum efficiency 
with time and affects all wavelengths, particularly those 
above 6 A. The Chandra instrument teams have studied 
the spectral profile of the contaminant absorption us- 
ing LETGS spect ra and have modeled its structure and 
time-dependence l|Marshall et al.ll2004|) . An estimated 
correction for the contaminant is now available through 
the Chandra website 3 . The biggest effect is an overall 
reduction in effective area due to the large optical depth 
of carbon. In addition, small edges from oxygen and flu- 
orine are also detected. The ARFs used in this analysis 
have been corrected for the effect of the contaminant. 

The plus and minus first order spectra and response 
files were summed to provide a single first order spectrum 
and response. Background spectra were extracted from 
the standard LETG background regions. After back- 
ground subtraction, the data were binned to a resolution 
of 0.075 A to ensure good statistics. Spectral analysis of 
the Chandr a observation of 4U 1850 —087 was performed 
using ISIS l|Houck fc Denicolall200C|) . 

2.2. Chandra Observation of 4U 1822- 000 

Chandra observed 4U 1822-000 on 2003 August 3 
for 2 ks using the High Energy Tra nsmission Grat- 
ing S pectrometer (HETGS) and ACIS l|Weisskonf et alJ 
|2003) . The HETGS carries two transmission gratings: 
the Medium Energy Gratings (MEGs) with a range of 
2.5-31 A (0.4-5.0 keV) and the High Energy Gratings 
(HEGs) with a range of 1.2-15 A (0.8-10.0 keV). The 
HETGS spectra are imaged by ACIS, which provides 
both an undispersed (zeroth order) image and dispersed 
spectra from the gratings. The various orders overlap 
and are sorted using the intrinsic energy resolution of 
the ACIS CCDs. The first-order MEG (HEG) spectrum 
has a spectral resolution of AA = 0.023 A or i?=110-1350 

2 http://asc.harvard.edu/ci ao/ 

3 See |http://asc. harvard.edu/ciao/threads/a ciscontam/ 
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(0.012 A or #=100-1250). The detected zeroth and first 
order count rates for 4U 1822-000 were 0.7 and 11.6 
cts s _1 , respectively. 

The level 1 event file for the 4U 1822-000 observa- 
tion was processed using CIAO v3.0. The zeroth or- 
der Chandra spectrum of 4U 1822—000 was heavily af- 
fected by pileup (pileup fraction >80%) and was not 
used in this analysis. We checked the dispersed spec- 
tra of 4U 1822—000 and found that pileup was minimal 
(<10%). The standard CIAO tools were used to cre- 
ate ARFs and RMFs for the MEG and HEG +1 and 
— 1 order spectra. The ARFs include the contamination 
correction mentioned previously. The ARFs were com- 
bined when the plus and minus first order spectra were 
added for the MEG and HEG separately. We created 
background spectra for the MEG and HEG using the 
standard background regions. The data were binned to 
ensure at least 50 counts per bin. The Chandra spectral 
analys is of 4U 1822 -000 was performed using XSPEC 
vll.3 ijArnau 
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2.3. XMM-Newton Observation of 4U 0513-40 

The XMM-Newton Observatory observed 4U 0513 
on 2003 April 1 for 24 ks. XMM-Newton carries 
three different instru ments, the European Photon Imag- 
ing C ameras (EPIC; iStriider et all 120011: iTurner et all 
20011). the Reflection Grating Spectrometers (RGS; 



den Herder et aUl200lD . and the Optical Monitor (OM; 
Mason et al.ll200lT) . The OM data was not used in this 



analysis. The EPIC instruments consist of three CCD 
cameras, MOS 1, MOS 2, and pn, which provide imag- 
ing, spectral, and timing data. The pn camera was run 
in small window mode and the two MOS cameras were 
run in full window mode. The EPIC instruments pro- 
vide good spectral resolution (A£;=50-200 eV FWHM) 
over the 0.3-12.0 keV range. There are two RGS on- 
board XMM-Newton which provide high-resolution spec- 
tra (i?=100-500 FWHM) over the 5-38 A (0.33-2.5 keV) 
range. The grating spectra are imaged onto CCD cam- 
eras similar to the EPIC-MOS cameras which allows for 
order sorting of the high-resolution spectra. The de- 
tected pn count rate for 4U 0513—40 was 20.8 cts s" 1 , 
while the first order combined RGS count rate was 2.3 
cts s _1 . 

The XMM-Newton data were reduced using the Science 
Analysis System (SAS) version 5.4.1. Standard filters 
were applied to all XMM-Newton data. The EPIC-pn 
data were reduced using epchain. The pn spectrum was 
extracted from a region 45" in radius around the source. 
Response files were generated for the pn spectrum us- 
ing rmfgen and arfgen. We reduced the EPIC-MOS 
data using emchain. The MOS data were found to have 
considerable pileup and were therefore not used in this 
analysis. The RGS data were reduced using rgsproc, 
which produced standard first order source and back- 
ground spectra and response files for both RGS. We 
grouped the pn spectrum to oversample the energy res- 
olution of the CCD by no more than a factor of three. 
The RGS spectra were grouped to ensure that each bin 
had at least 20 counts. All XMM-Newton spectr al analy- 
sis of 4U 0513—40 was performed using XSPEC iJArnaudl 
1996). 

3. ANALYSIS AND RESULTS 




14 16 18 20 22 24 
Wavelength (A) 

Fig. 1. — (Top panel) Combined LETG first order spectrum 
of 4U 1850—087 with best-fit model including the oxygen, iron, 
and neon edge models. (Bottom panel) Fractional residuals ([data- 
model] /model) of the spectral fits shown above. The absorption 
lines at 23.5 and 23.3 A are the expected interstellar absorption 
lines from O I and O II. There are no other features detected 
which are consistent with narrow emission or absorption lines. The 
arrows mark the neon, iron, and oxygen absorption edges at 14.3, 
17.5, and 23 A respectively. 



3.1. Chandra Observation of 4U 1850-087 

We previously showed that 4U 1850—087 was one of 
four sources to have a feature at 0.7 keV in its ASCA 
spectrum which could be explained by an unusual Ne/O 
ratio in the absor bing material along the line of sight 
ijJuett et alJl200if) . High-resolution spectroscopy of the 
other three sources has confirmed the unusual Ne/O 
ratios l|Paerels et all 120011: Ouett fc Chakrabartvll2003|) . 
To investigate the absorption along the line of sight to 
4U 1850-087, we fit the first order LETGS spectrum 
over the 1.5-25.5 A wavelength range with a model that 
separately fit for the column densities of the oxygen, iron, 
and neon edges. The 5-10 A wavelength range was ex- 
cluded from the fits in order to avoid the piled up region 
of the spectrum. All errors are 90%-confidence unless 
otherwise noted. 

The continuum model consisted of a power law + 
blackbody and absorption by the tbvarabs model with 
the abundances for oxygen, iron, and neon set to zero 
and all othe r s equ al to the interstellar abundances of 
iWilms et alJ ((2000) . The equivalent hydrogen column 
density Ah in tbvarabs was set to the expected in- 



reddening £b-v to the cluster l)Kuulkers et alJ Lww.^ 
and using the rela tionship between Eb-v and Ah of 
iPredehl fc SchmittJ (|1995|) . The AT-shell edges of oxy- 
gen and neon were fit using a standard edge model, 
while the Fe-L edge was fit using a cust om multiplicative 
model based on the cross-sections of iKortright fc Kiml 
(2000). Two Gaussian absorption lines were included 
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TABLE 1 

Best-Fit Continuum Spectral Parameters 51 



Source 




r 


AS 


kT bb (keV) 


^-km/^lO kpc 


Flux c 


f d 
/pi 




4U 1850- 
4U 0513- 
4U 1822- 


087 

40 

000 


1.48±0.05 
1.963±0.017 
1.3±0.3 


1.32±0.09 
2.04±0.06 
8±2 


0.479±0.009 
0.75±0.08 
0.86^ 


113±13 
2.0±0.7 
4&±? 7 


1.45 
0.97 
8.4 


0.66 
0.96 
0.73 


1.3/242 
1.12/1612 
1.47/452 



a All errors quoted at the 90%-confidence level. 

b Power-law normalization in units of 10~ 2 photons cm -2 s — 1 keV - 
c The absorbed 0.5-10 keV flux in units of 10 — 10 erg cm -2 s — 1 . 
d Fraction of the absorbed flux in the power-law component. 



at the oxygen edge to represent the ls-2p absorption 
lines from O I and O II. For oxygen and neon, the 
atomic absorption cross -sections were taken from the the- 
oreti cal calculations of iGorczvca fc McLaughlin! ( 2000) 
and iGorczvcal l|2000j) . respectively. The spectrum of 
4U 1850—087 and best-fit spectral model including the 
absorption edges is shown in Figure ^ The best-fit 
power-law and blackbody parameters are given in Ta- 
ble □ The absorbed 0.5-10 keV flux for 4U 1850-087 
was 1.45 x 10~ 10 erg cm~ 2 s _1 with the power-law 
component contri buting 66% of the flux. At a dis- 
tance of 8.2 kpc ijKuulkers et alJfeOOSj) . the absorption 
corrected 0.5-10 keV luminosity of 4U 1850—087 was 
1.35 x 10 36 erg s , assuming isotropic emission. 

The best-fit absorption edge and line values are given 
in Tabled The best-fit edge positions are consistent with 
the expected positions of the edges. We find a Ne/O ra- 
tio along the line of sight to 4U 1850-087 of 0.22±0.05, 
consistent with the expect ed interstellar value of 0.18 us- 
ing the ISM abundances of lWilms et al.1 fl200flfl. as well as 
that found in the XMM-Newton spectrum (iSidoli et al.l 
120041) . The oxygen column density implies a hydrogen 
column density of Ah = (4.2 ± 0.4) x 10 21 cm" 2 as- 
suming standard ISM abundances l)Wilms et alJl2000T) in 
agreement with the predicted Ah from the neon and iron 
edges, although a factor of two larger than the expected 
Ah of 2 x 10 21 cm" 2 fro m the reddening to the clus- 
ter ijKuulkers et al.lf2003[l . Other than the interstellar 
oxygen absorption lines, no other emission or absorption 
lines were apparent in the spectrum of 4U 1850—087. We 
performed a careful search of the Chandra spectral resid- 
uals to place limits on the presence of any spectral fea- 
tures in the 12-25 A wavelength range. Gaussian models 
with fixed FWHM = 2000 km s" 1 were fit at each point. 
(This width is similar to that see n in emission lines from 
other LMXBs, e.g., 4U 1626-67. ISchulz et al.ll200ll and 
EXO 0748-67. ICottam et al.ll2001ft . Frorn~this. we can 
place a 3c upper limit on the E W of any line feature, 
either emission or absorption (see iProtassov et al.ll2002l 
for a cautionary note). The EW limit increases with 
wavelength, varying from 0.05 A at 12 A to 0.25 A at 
24 A. 

The LETGS first order dispersed spectrum of 
4U 1850—087 has an average count rate of 2.80 ± 
0.17 cts s _1 . We examined the total count rate, as well 
as the count rates in two different energy ranges, to check 
for changes in the spectral state. We found no evidence 
for any change of state during the Chandra observation. 
To look for periodic modulations of the X-ray flux, we 
created lightcurves from event files after barycentering 



and randomizing the event arrival times. (Randomiz- 
ing of the event arrival times consists of adding a ran- 
dom quantity uniformly distributed between zero and the 
readout time of 0.74 s in order to avoid aliasing caused by 
the readout time.) We searched for modulations of the 
X-ray flux with frequencies between lxlO -5 and 5xl0 -2 
Hz. We found no evidence for periodic modulation with 
a 90%-confidence upper limit of 1.3% for the fractional 
rms amplitude. 

To determine the position of 4U 1850—087, we first 
applied the aspect offset correction available at the 
Chandra website 4 . Using the CIAO tool celldetect, 
the zeroth order source position for 4U 1850—087 was 
determined: R.A.=18 h 53 m 04!86 and Dec=-8°42'20"4, 
equinox J2000.0 (90% confidence error radius of 076). 
This position is Of! 6 from the optic al counterpart posi- 
tion given bv I Anderson et all |jTp93). 

3.2. XMM-Newton Observation of 4U 0513-40 

To determine the appropriate continuum model for 
4U 0513—40, we fit the pn and RGS spectra simultane- 
ously with either an absorbed power-law or an absorbed 
power-law + blackbody model. The pn spectrum was 
fit over the energy range 1.5-12 keV, while the two RGS 
spectra were fit over the range 0.34-2.0 keV. Residuals on 
the order of 10% were seen in the pn at energies between 
0.5-1.5 keV. While this is consistent with estimates of 
the effective area calibration for the pn, given the high 
count rate of 4U 0513—40, these residuals resulted in in- 
flated xt values so we chose to exclude this region of 
the pn spectrum from our fits. A constant was also 
included to allow for normalization differences between 
the pn and RGS. The absorbed power-law + blackbody 
model was a significantly better fit (significance of 98% 
as calculated by an F-test) to the data and we adopt 
this as the appropriate continuum model. The best-fit 
power-law and blackbody parameters are given in Ta- 
ble □ The absorbed 0.5-10 keV flux for 4U 0513-40 
was 9.66 x 10 -11 erg cm -2 s _1 with the power-law com- 
ponent co ntributing 96% of th e flux. At a distance of 
12.1 kpc ijKuulkers et alJ l2003). this gives a luminosity 
of 1.8 x 10 36 erg s" 1 . 

To determine the column density along the line of sight 
to 4U 0513—40, we fit the RGS spectra with separate 
models for the absorption edges from oxygen, iron, and 
neon as was done for 4U 1850—087. The continuum 
model was fixed to that found in the pn and RGS com- 
bined fit. A Gaussian absorption line was included to 

4 See |http://asc.harvard.edu/cal/ASPECT/celrjaon/index.html| 
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Fig. 2. — (Top panel) Combined RGS first order spectrum of 
4U 0513—40 with best-fit model including the oxygen, iron, and 
neon edge models. (Bottom panel) Fractional residuals ([data- 
model] /model) of the spectral fits shown above. The absorption 
line at 23.5 A is the expected interstellar absorption line from O I. 
There are no other features detected which are consistent with nar- 
row emission or absorption lines. The arrows mark the neon and 
oxygen absorption edges at 14.3 and 23 A respectively. 



model the ls-2p absorption line from O I. The spec- 
trum of 4U 0513—40 and best-fit spectral model are 
shown in Figure The best-fit absorption line and 
edge parameters are given in Table [3 The absorption 



to 4U 0513-40 is very low (N H 



10 



20 



cm 



') mak- 



ing a determination of both the wavelength and depth 
of the iron and neon edges difficult. We therefore fixed 
the iron and neon edge wavelengths to the expected val- 
ues of 17.5 and 14.3 A, respectively. Given the strength 
of the oxygen edge complex, its position was allowed 
to vary and the best-fit edge wavele ngth is comp arable 
to what is found in other LMXBs il.liiott el all (2004). 
The oxygen edge measurement implies an equivalent hy- 
drogen column of TVh = (7.6 ± 1. 2) x 10 20 cm" 2 as - 
suming standard ISM abundances (j Wilms et al .1 12000). 
The neon and iron edges are consistent with this value. 
The X-ray derived TVh value is greater than the value of 
TVh = 1 .1 x 10 20 cm" 2 derive d from the reddening to the 
cluster l)Kmilkers et al.l2 003). No emission or absorption 
lines were apparent in the spectrum of 4U 0513—40. We 
performed a search of the XMM-Newton RGS spectral 
residuals over the 6-36 A wavelength range. The 3a up- 
per limit on the EW of any spectral feature is roughly 
constant at 0.1 A over the range 6-30 A and then in- 
creases to 0.4 A at 35 A for Gaussian models with fixed 
FWHM = 2000 km s" 1 (see iProtassov et al.ll2002l for a 
cautionary note). 

The 4U 0513—40 observation had average count rates 
of 20.8±1.2 and 7.1±0.6 cts s" 1 for the pn and com- 
bined RGS, respectively. Using the pn data, we checked 
for changes in the spectral state by examining the total 



count rate, as well as the count rates in two different en- 
ergy ranges. There was no evidence for a change in the 
spectral state of 4U 0513-40 during the XMM-Newton 
observation. Due to the low count rate of this obser- 
vation (-C 1 count per readout), the low frequency end 
of t he power spectrum is affected by instrumental noise 
(see lKuster et alJ 120021 and references therein). To re- 
move the noise, we calculated the average power at each 
frequency using an 11 bin moving average and then di- 
vided by this value. We verified that the resulting pow- 
ers were consistent with an exponential distribution, and 
then searched for modulation of the X-ray flux with fre- 
quencies between 3 x 10~ 4 and 1 Hz. No modulation 
was detected with a 90%-confidence upper limit on the 
fractional rms amplitude of 0.7%. 

3.3. Chandra Observation of 4U 1822- 000 

We fit the Chandra HEG and MEG first-order spec- 
tra jointly with both an absorbed power-law and an ab- 
sorbed power-law + blackbody model over the wave- 
length ranges 1.8-11.3 A for the HEG and 2.1-14.3 A 
for the MEG. The data were binned to have at least 50 
counts per bin which given the low number of counts 
from the source resulted in significantly reduced spectral 
resolution. This combined with the high Nn along the 
line of sight resulted in few counts being detected short- 
ward of 0.9 keV. Therefore, it was not necessary to use 
the complex absorption model presented ear lier. Instead, 
we used the in terstellar absorption model of lWilms et alJ 
(tbabs; 2000) with the abundances in XSPEC set to 
wilm and the cross-section set to vern. The power-law 
+ blackbody model produced a significantly better chi- 
squared value than the power-law alone (significance of 
98.4% as calculated by an F-test), so we take a power- law 
+ blackbody model as the appropriate continuum model 
for the Chandra spectra. The best-fit equivalent hydro- 
gen column density was iVn = (0.97 ± 0.18) x 10 22 cm~ 2 
which compares well with the value derived from infrare d 
dust maps (iV H ~ 0.8 x 10 22 cm- 2 : ISchlegel et alJll99^ . 
The best-fit power-law and blackbody parameters are 
given in Table Q The absorbed 0.5-10 keV flux for 
4U 1822-000 was 8.4 x lO" 10 erg cm" 2 s" 1 with the 
power-law component contributing 73% of the flux. 

The combined MEG + HEG first order dispersed spec- 
trum of 4U 1822-000 has an average count rate of 
11.3 ± 0.5 cts s -1 . We examined the total count rate, 
as well as the count rates in two different energy ranges, 
to check for changes in the spectral state. We found 
no evidence for any change of state during the Chan- 
dra observation. To look for periodic modulations of the 
X-ray flux, we created lightcurves from event files after 
barycentering and randomizing the event arrival times. 
To increase statistics, we included events from the first 
through tenth grating orders. We searched for modula- 
tions of the X-ray flux with frequencies between 2xl0 -4 
and 5x 10~ 2 Hz. We found no evidence for periodic mod- 
ulation with a 90%-confidence upper limit of 2.6% for the 
fractional rms amplitude. 

To determine the position of 4U 1822—000, we first 
applied the aspect offset correction. Using the CIAO 
tool tgdetect, the zeroth order source position for 
4U 1822-000 was determined: R.A.=18 h 25 m 22!02 and 
Dec=-00 o 00'43"0, equinox J2000.0 (90% confidence er- 
ror radius of 076). This position is 0'.'9 away from the 
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TABLE 2 

Best-Fit Absorption Feature Parameters 51 



Feature 


A (A) N x (10 18 cm" 2 ) 


EW (mA) 




4U 1850- 


-087 




O-K edge 
O I ls-2p line 
O II ls-2p line 
Fe-L edge 
Ne-K edge 


22.79+0"! 

23.511 ±0.004 

23 363+ 04 
zo.ouo_ 014 

17.52 ±0.07 

14.40 ±0.08 


2.1 ±0.2 

0.08 ±0.03 
0.46 ±0.10 


70 ±30 
40 ±40 




4U 0513- 


-40 




O-K edge 
O I ls-2p line 
Fe-L edge 
Ne-K edge 


9^ n r >+ u-l; - i 

zo.uo_ gg 
23.4±g;| 
17.5 (fixed) 
14.3 (fixed) 


0.37 ±0.06 

< 0.014 
0.10 ±0.08 


31 ± 26 



a All errors quoted at the 90%-confidence level. 



optical counterpart position ijGottwald et alJ ll991). 
3.4. Chandra Observation of 4U 1905+000 

We observed 4U 1905+000 on 2003 March 4 for 5 ks 
with Chandra in the HETG/ACIS configuration. We did 
not detect a source during this observation. From the 
non-detection, we place an upper limit on the flux of 
1 x 10~ 14 erg cm" 2 s" 1 assuming a 0.3 keV blackbody 
and Nf\ = 1.9 x 10 21 cm" 2 . From a radius expansion 
burst. IChevalier fc Ilovaiskvl l) 19901) derived a distance to 
4U 1905+000 of 8+1 kpc. At this distance, the upper 
limit on the source luminosity is 1 x 10 32 erg s _1 , making 
4U 1905+000 one of the faintest known quiescent NS 
systems. 

In an effort to determine when the source went into 
quiescence, we retrieved the archival data from ROSAT 
and ASCA observations of 4U 1905+000. ROSAT obser- 
vations of Aql X-l would have included 4U 1905+000 
in the field of view. We retrieved the two PSPCB 
observations of Aql X-l performed on 1992 October 
15 and 1993 March 24. No source at the position of 
4U 1905+000 was found with a upper limit on the flux 
of 6 x 10~ 13 erg cm~ 2 s _1 . Similarly, the ASCA ob- 
servation of 4U 1905+000 performed on 1995 April 09 
found no source at the position of 4U 1905+000. From 
this work, we conclude that 4U 1905+000 went into its 
quiescent state in the late 1980s or early 1990s after an 
extended outburst phase. 

4. DISCUSSION 

In this paper, we have shown that the Chan- 
dra/LETGS spectrum of 4U 1850-087, the XMM- 
Newton/RGS spectrum of 4U 0513-40, and the Chan- 
dra/HETGS spectrum of 4U 1822-000 are all well fit 
by an absorbed power-law + blackbody model. No un- 
usual abundance ratios are implied by the detected ab- 
sorption features, and no emission or absorption features 
were found other than the expected absorption lines due 
to interstellar oxygen. For 4U 1850—087, this is sig- 
nificantly different from what is implied by the earlier 
ASCA observation, which had established 4U 1850-087 
as being part of a class of four NS/LMXBs all having a 
similar feature at 0.7 keV in their low-resolution spectra 
ij.Iuett et alJl200ll) . This feature had originally been at- 
tributed to u nresolved line emission from iron and oxy - 
gen (see e.g., iChristian etaD 11994 iWhite et al]ll997|) . 



However, a high-resolution observation of the brightest of 
these sources, 4U 0614+091, with Chandra, failed to de- 
tect any emission lines, finding instead an unusually high 
Ne/O number ratio in th e absorption along the line of 
sight l|Paerels et al.ll200lj) . We had subsequently pointed 
out that the ASCA spectra of all four sources are well fit 
without a 0.7 keV emission line using a model that in- 
cludes photoelectric absorption due to excess neon along 
the lines of sight and presumably local to the sources 
H-Tuett et al.l l200T|). 

Of the other four sources that we previously pro- 
posed as neon-rich systems based on their low-resolution 
ASCA spectra, three (4U 0614+091, 2S 0918-549, 
and 4U 1543—624) have also shown unusual Ne/O 
abundance ratios in their high-resolution X-ray spectra 
(|Paerels et al.ll2001l ITuett fc Chakrabartvll2003[) . How- 
ever, a comparison of the Chandra, XMM-Newton, and 
ASCA results for 2S 0918-549 and 4U 1543-624 also 
revealed variability in the measured Ne/O ratio, which 
seemed to be assoc iated with changes in the so urce con- 
tinuum properties ijJuett fc C hakrabartv 2003). We pre- 
viously suggested that changes in the spectral proper- 
ties of the sources could lead to ionization effects that 
would produce the variable Ne/O ratio found for the 
neutral edges. Ionization would affect oxygen more than 
neon, inflating the measured neutral Ne/O ratio. Un- 
fortunately, this means that the measured Ne/O ratio 
would not reflect the donor composition. 

Here we propose that for 4U 1850—087 as well, ion- 
ization effects due to continuum spectral variations are 
also responsible for the differences between the ASCA 
and Chandra results. When the Chandra spectrum is 
rebinned to a resolution comparable to that of ASCA, 
the Chandra data do not show the 0.7 keV feature seen 
previously with ASCA (see Figure [3J|. This clearly in- 
dicates that the source spectrum changed between 1995 
and 2002. Given our interpretation of the feature as be- 
ing due to an unusual Ne/O ratio in the absorbing ma- 
terial, the lack of the 0.7 keV feature in the Chandra 
spectrum is consistent with the measured Ne/O ratio of 
0.22+0.05, which is close to the interstellar value. Had 
the feature still been present in the binned Chandra spec- 
trum, then our interpretation would have been incorrect. 
The variation of the Ne/O ratio in the absorbing material 
is accompanied by differences in the continuum spectral 
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Fig. 3. — Comparison of the ASCA (top panel) and Chandra 
(bottom panel) spectrum for 4U 1850—087 plotted at the same 
energy resolution. Shown are the fractional residuals ([data — 
model] /model) of an absorbed power law + blackbody model. The 
ASCA spectrum shows a residual feature at 0.7 keV (marked by an 
arrow) which is absent in the Chandra data. This feature was pre- 
viously attributed to an unusual Ne/O abundance ratio in the ab- 
sorption. The lack of the 0.7 keV feature in the Chandra spectrum 
supports our suggestion that the source spectrum has changed be- 
tween the two observations. The residual in the Chandra spectrum 
at 0.55 keV is due to misalignment of the oxygen edge between the 
data and model. 



properties of 4U 1850—087. The Chandra observation of 
4U 1850—087 was taken during a lower luminosity state 
than the ASCA observation which showed the unusual 
Ne/O ratio. This is comparable to what was observed in 
the other proposed neon-rich sources 2S 0918—549 and 
4U 1543-624. A more detailed analysis of the XMM- 
Newton data may shed more light on this source. Given 



that the XMM-Newton results sho w a Ne/O ratio sim ilar 
to that found in the Chandra data ijSidoh et all2004D . we 
would expect that the continuum properties should also 
be similar. 

While X-ray spectra may be inconclusive in determin- 
ing the donor composition of these systems, new results 
from optical spectroscopy have confirmed their unusual 
nature. Optical spectra of three candidate ultracompact 
binaries with suggested neon-rich do nors. 4U 0614+091 . 
4U 1543-624, and 2S 0918-549 (I.Tuett et all l200l|) . 
showed no prominent hydrogen or h elium emission lines 
l|Nelemans et all 120041 IWand 120041) . In the b rightest 
of these, 4U 0614+09l HNelemans et~aTl l|2004l) found 
emission lines which they attributed to carbon and oxy- 
gen. The other spectra were not as high quality, but 
hints of similar spectral features were seen. The op- 
tical spectra seem to confirm our su ggestion that the 
donors in these system s are C-0 WDs (I.Tuett et al.ll20rrl 
INelemans et alJ I2004|) . The notion of a class of ul- 
tracompact LMXBs with neon-rich C-0 WD donors is 
motivated by th e definite detection of one such sys - 
tem, 4U 1626-67 (|Schulz et alJ2001tlHomer et al.l 2002). 
However, not all ultracompact binaries are expected to 
have C-0 donors. The burst properties of the 11-min 
binary 4U 1820—30 support a helium WD donor for that 
system, while the properties of the 50-min X-ray dip- 
per 4U 1916—05 point to a hydrogen-deficient but not 
yet degenerate donor. As more observations allow for 
the donor compositions of the ultracompact binaries to 
be determined, the results will provide a unique test of 
binary evolution theory. Multiwavelength observations 
may be the best way of determining the donor proper- 
ties in ultracompact binary systems. 
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